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ABSTRACT 


A  group  of  stations  in  the.  North.  American  Arctic  region  have 
been  analyzed  for  statistical  determination  of  temperatures  at 
mandatory -pressure  levels.   For  each  station  the  temperature  at  a 
key  level,  peculiar  to  that  station,  has  been  forced  in  at  the 
first  step  and  retained  at  each  subsequent  step,  in  the  development 
of  stepwise  regression  equations  which  predict  temperatures  at  the 
mandatory  levels.   In  general,  eight-step  predictions,  in  terms  of 
inter-level  thicknesses,  were  found  to  give  optimum  specification 
of  the  desired  temperatures.   The  best  estimate  of  the  regional 
atmosphere,  which  is  conditionally  dependent  upon  the  existence  of 
an  extreme  1%  probability  of  the  forcing-level  temperature,  is 
obtained  with  a  high  degree  of  confidence. 
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I.   INTRODUCTION 

Although  many  factors  must  be  considered  in  determining  the 
thermal  environment  of  airborne  ordnance,  the  most  dominant  are  the 
flight  conditions  and  the  atmospheric  temperature  profile.   While  air- 
craft and  missile  performance  capabilities  are  well  defined,  there  is 
a  great  deal  of  variability  in  the  atmospheric  temperature  profile 
which  may  be  traversed  by  the  vehicle.   Therefore  many  model  atmo- 
spheres have  been  developed  for  meteorological  purposes,  but  none  is 
completely  satisfactory  for  use  in  aerothermodynamic  heating 
calculations. 

Supplemental  standard  atmospheres  have  been  developed  by 
COESA(1966)  to  provide  information  on  spatial  and  temporal  variability 
across  selected  latitude  bands  and  by  seasons.   However  the  supplement 
standards  are  intended  to  be  means  with  respect  to  the  latitude  bands 
and  seasons,'  and  are  not  representative  of  the  extremes  which  must  be 
considered  for  aerodynamic  design  purposes.   Richard  and  Snelling  (1971) 
proposed  the  adoption  of  a  revised  MIL-STD-210  hot  and  cold  model  atmo- 
spheres which  represent  global  extreme  values  at  successive  2  km  eleva- 
tions.  These  extremes  are  deviations  from  the  station  mean  at  the 
elevation  involved,  which  are  exceeded  at  the  1%,  5%  and  10%  probability 
levels,  and  were  derived  through  analysis  of  extreme  data  from  global 
meteorological  temperature  maps  analyzed  at  constant  pressure  levels. 
More  recently  a  second  draft,  MIL-STD-210B,  has  been  published  (1972). 
Both  of  these  extreme  sets  of  profiles  pose  serious  design  problems  when 
changes  of  elevation  are  involved. 


The  most  important  limitation  of  the  two  recently  proposed 
extreme-atmosphere  sets  is  that  they  cannot  be  realistically  utilized 
in  conjunction  with  flight  profiles  and  trajectories  which  involve 
changes  in  altitude.   The  temperature  extremes  at  differing  altitudes 
are  in  most  cases  representative  of  different  geographical  locations 
rather  than  of  conditions  occurring  in  real  time  over  a  real  location. 
Polar  and  tropical  standard  and  extreme  atmospheres  have  been  developed 
for  such  applications  but  are  still  subject  to  the  deficiencies  noted 
above. 

TABLE  1.   Proposed  locations  of  the  1%  world-wide  extreme  temperatures 
at  indicated  levels  (AFCRL,  1970), 


Pressure 

Level 

altitude 

1%  Cold  extreme 

1%  Warm  extreme 

(ft) 

location 

location 

SFC 

0 

Oymyakon,  USSR 

Insalah,  Algeria 

850  mb 

4780 

Oymyakon,  USSR 

Insalah,  Algeria 

700  mb 

9882 

Hall  Beach,  NWT  » 

Babylon,  Iraq 

500  mb 

18287 

Resolute,  NWT    * 

New  Delhi,  India 

300  mb 

30066 

Thule,  Greenland* 

New  Delhi,  India 

200  mb 

38661 

Thule,  Greenland" 

Alert,  NWT      * 

150  mb 

44646 

Karachi,  Pakistan 

Alert,  NWT      * 

100  mb 

53084 

Singapore 

Thule,  Greenland* 

The  object  of  this  study  was  to  develop  the  most  probable  verti- 
cal profile,  conditional  with  the  existence  of  a  1%  extreme  at  a  given 
level,  for  each  station  set  forth  in  Table  1  (but  here  limited  to  those 
marked  by  an  asterisk).   Each  station  of  Table  1  has  been  associated 
with  a  listed  1%  extreme  for  the  indicated  month  according  to  a  pre- 
liminary draft  of  the  MIL-STD-210B  proposed  atmosphere,  which  is 
graphically  illustrated  in  Fig.  1.   The  asterisk-marked  stations  nil 
occur  :i  n  the  Canadian  Arctic  and  are  either  extremely  cold  at  the 
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indicated  level  in  winter  or  extremely  warm  in  the  summer  stratosphere. 
This  subset  of  the  stations  in  Table  1  was  isolated  from  the  full  set 
because  the  data  were  readily  available  in  conveniently  summarized 
form  from  checked  rawinsonde  data.   The  financial  assistance  and  kind 
cooperation  of  the  Commander,  Naval  Weather  Service  is  gratefully 
acknowledged  in  acquiring  the  data  on  magnetic  tape  from  the  National 
Environmental  Data  Service,  Ashville,  North  Carolina.   For  each  station 
and  month  of  interest,  four  years  (1967-70)  of  atmospheric  soundings 
were  processed  to  determine  the  temperature  means,  standard  deviations 
and  inter- level  correlations  required  for  application  of  the  multiple 
regression  analysis  procedure  to  be  used  in  the  determination  of  the 
most  probable  temperature  profile  at  the  1%  probability  level. 
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II.   DATA  PROCESSING 

For  each  of  the  stations  marked  with  an  asterisk  in  Table  1,  four 
years  of  checked  rawinsonde  data,  for  the  period  1967  through  1970, 
was  provided  by  the  Naval  Weather  Service  Environmental  Detachment, 
Ashville.   Those  stations  listed  under  the  cold  extreme  of  Table  1 
were  subject  only  to  January  data  samples,  and  those  under  the  warm 
extreme,  only  to  July  samples.   The  individual  sounding  samples  were 
arranged  in  the  sequential  format  indicated  in  Table  2  and  recorded 
on  magnetic  tape.   Height  and  temperature  data  for  the  indicated 
pressure  levels  were  arranged  with  six  sets  in  each  of  the  first  three 
rows  and  five  sets  in  the  fourth  row.   Note  that  the  data  is  arranged 
at  50  mb  increments  from  1000  to  200  mb  inclusively  and  at  somewhat 
smaller  increments  thereafter. 


TABLE  2.   Arrangement  of  sounding  temperature  and  geopotential 
height  data  by  pressure  levels  k=l,...,23  in  each  sounding. 


Index 
level 

k  = 
(mb) 

1 
1000 

2 
950 

3 
900 

4 
850 

5 
800 

6 
750 

Index 
level 

k  = 
(mb) 

7 
700 

8 
650 

9 
600 

10 
550 

11 
500 

12 
450 

Index 
level 

k  = 
(mb) 

13 
400 

14 
350 

15 
300 

16 
250 

17 
200 

18 
175 

Index 
level 

k  = 
(mb) 

19 
150 

20 
125 

21 
100 

22 
80 

23 

70 

With  the  data  in  the  format  of  Tabic  2,  it  is  convenient  to 
determine  the  mean  temperature  T(P_)  and  standard  deviation  a      for  the 
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indicated  station  levels  specified  by  Table  1.   The  temperature  T(P  ) 
is  to  be  used  in  a  multiple  regression  analysis  to  determine  the  most 
probable  vertical  temperature  profile  at  all  mandatory-sounding  levels 
associated  with  that  "forcing-level  temperature"  T(P  ).   A  subset  of 
such  station  data  is  further  analyzed  for  determination  of  the  above 
statistical  parameters  for  those  soundings  having  10%  and  1%  extreme 
temperatures  at  the  forcing  level.   The  T  (.10)  and  T  (.01)  extreme 
temperatures  can  then  be  forced  into  the  derived  full-sample  regression 
procedure  to  obtain  the  corresponding  extreme  profiles. 

TABLE  3.   Nominal  temperature  extremes  at  indicated  stations  and 
pressure  levels  at  the  10%  and  1%  values  of  a  normal  distribution. 


January 

Cases 

Full 

Sample 

Forcing 

Sample 

Size  10% 

Station 

Level  T 

Tj(.lO') 

TjC.01) 

Size 

Extreme 

Hall  Beach 

T(700) 

-35.6 

-42.3 

243 

26 

Resolute 

T(500) 

-47.1 

-51.3 

237 

25 

Thule 

T(300) 

-62.0 

-67.3 

213 

13 

Thule 

T(200) 

-63.5 

-70.6 

213 

18 

July 

Cases 

Thule 

T(100) 

-41.3 

-39.8 

236 

23 

Alert 

T(150) 

-39.7 

-37.9 

244 

22 

Alert 

T(200) 

-38.6 

-35.7 

244 

18 

Table  3  lists  the  forcing  level  pressures  and  sample  sizes  with 
the  10%  and  1%  limiting  temperatures  for  each  station  under  considera- 
tion.  Estimates  of  the  temperatures  T(P  ),  at  the  probability  extri 
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of  10%  and  1%  have  been  computed  assuming  a  Gaussian  distribution  at 
each  J-level.   These  estimates,  T  (.10)  and  T  (.01)  of  Table  3,  are 
derived  using  the  computed  full  sample-means  T (P  )  and  standard  devia- 
tions a   of  the  1967-70  samples  together  with  the  well-known 
relationships 

T  (.10)  =  T  -  1.2817a 

_   ,  (1) 

T  (.01)  =  T  '  -  2.3267CJ 

where  the  plus  sign  is  to  be  used  for  a  warm  extreme  and  the  minus  sign 

for  a  cold  extreme. 

Analysis  of  the  data  revealed  an  insufficient  number  of  soundings 

occurring  at  the  1%  extreme  of  T(P  )  for  a  test  of  the  statistical 

significance  of  these  results.   It  was  decided  therefore,  to  give 

primary  consideration  to  the  10%  extreme  T (P  )  soundings  and  to  deduce 

the  properties  of  the  1%  subset  as  a  limiting  case  of  the  10%  sample. 

Thus  the  same  multiple  regression  treatment  used  to  predict  the  T(P  ) 

mandatory-level  temperatures  for  the  full  sample  has  been  performed 

upon  the  "10%  subset"  containing  T (P  )  extremes. 


15 


III.   THE  REGRESSION  METHOD 

A.   THE  STEPWISE  REGRESSION  PROCEDURE 

Many  problems  in  research  require  the  extensive  analysis  of  large 
amounts  of  data,  and  the  process  for  handling  it  should  be  made  as 
automatic  and  rapid  as  possible.   To  this  end,  the  Biomedical  Computer 
Programs  (Dixon,  1966)  were  developed.   One  of  these  programs,  BIMED 
02R,  provides  stepwise  regression  analysis  capability  and  will  be 
utilized  here. 

BIMED  02R  computes,  in  a  stepwise  manner,  a  sequence  of  multiple 
linear  regression  equations.   That  single  variable  is  added  to  the 
equation  at  each  step  which  produces  the  greatest  reduction  in  the  pre- 
vious unexplained  variance.   This  variable  is  also  the  one  which  has 
the  highest  partial  correlation  with  the  dependent  variable  at  the 
particular  step  in  the  analysis  of  the  variance.   Equivalent ly  it  is 
the  variable  which  would  give  the  highest  F-statistic  upon  entry  at 
that  step.   The  F-statistic  upon  entry  F  ,  is  expressed  at  step  r  as 


(Dixon  1966) 


where 


F  (l,n-r-l)  =  %  (C.E.V.,r)  -  %  (C.E.V.,r-l) 
r  %  (U.E.V. ,r) 


%  (C.E.V.,r)  is  the  percent  cumulative  explained 
variance  at  step  r, 

%  (C.E.V.,r-l)  is  the  percent  cumulative  explained 
variance  at  step  r-1, 

%  (U.E.V. ,r)  is  the  percent  unexplained  variance 
remaining  at  step  r. 
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This  study  utilized  a  statistical  model  expressible  in  the  form 

TM  =  co  +  ciT(V  +  eA 

Here  T   in  (2)  is  the  temperature  predictand  at  each  of  the  eight 
mandatory  levels  chosen  alternately  as 


(2) 


M 


V  V  V  V  V  V  V  V  and  T9 


but  with  T  5^  T  (P  )  ,  since  the  forcing- level  temperature  is  used  as  the 
forced  predictor  for  each  of  the  eight  mandatory  level  T   specifications 
for  the  station  under  consideration.   Thus  for  each  of  the  stations  in 
Table  3,  eight  of  the  following  nine  dependent  variables  were  specified 
for  use  in  the  multiple  regression  procedure: 

the  temperature  at  1000  mb ,  T(1000) 

the  temperature  at  850  mb ,  T(850) 

the  temperature  at  700  mb ,  T(700) 

the  temperature  at  500  mb ,  T(500) 

the  temperature  at  300  mb ,  T(300) 

the  temperature  at  200  mb ,  T(200) 

the  temperature  at  150  mb ,  T(150) 

the  temperature  at  100  mb,  T(100) 

the  temperature  at  70  mb,  T(70). 
T(P  )  is  forced  into  the  specification  of  each  of  the  eight  possible 
predictands  in  the  manner  symbolized  by  Eq.  (2).   X  ,  called  a  struc- 
ture function,  represents  a  linear  combination  of  up  to  seven  inde- 
pendent variables,  comprised  of  layer  thickness-values  and  the  derived 
coefficients  which  meet  the  requirements  of  the  stepwise  least  squares 
technique.   The  thickness  variables  Z,  ,  were  generated  within  the  BIMED 
02R  program  by  utilization  of  the  so-called  transgencration  technique 
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which  takes  the  difference  of  successive  tabulated  heights  in  a 
sequential  manner  from  the  surface  to  70  mb.   This  procedure  produced 
twenty-two  independent  thickness  variables  which  are  directly  related 
to  the  relevant  mean  virtual  temperature  in  the  layer  (P,  ,  P    ) ,  by 
the  hypsometric  equation.   These  thickness  predictors  proved  to  be 
highly  useful  in  the  specification  of  T  . 

In  the  specification  of  the  dependent  variable  T  ,  a  statement  is 
permitted  as  to  the  maximum  number  of  independent  variables  to  be 
tested  for  admission  to  the  regression  equation.   In  the  tests  con- 
sidered here,  an  upper  limit  of  eight  independent  variables  was  found 
in  general  to  give  optimal  specification  of  each  of  the  dependent 
variables.   Still  another  convenient  feature  of  the  BIMED  02R  program 
is  that  the  forced  entry  into  the  regression  of  the  temperature  value 
being  tested  as  an  extreme,  the  so-called  forcing-level  temperature, 
is  allowed  even  though  its  contribution  to  the  explained  variance  may 
become  small  at  the  final  step  in  the  regression  analysis. 

Tables  4  through  10  list  the  variables  selected  and  coefficients 
of  Eq.  (2)  for  each  station  extreme  level  considered.   The  listed 
structure  functions  X^  are  often  highly  dependent  upon  the  thickness 
variables  nearly  adjacent  to  the  level  represented  by  T  .   However, 
diversity  of  thickness-variable  selection  does  occur. 
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TABLE  4.   Full-sample  regression  equations,  with  associated  thickness- 
variables  at  Hall  Beach,  NWT,  for  the  specification  of  the  January 
temperatures  at  indicated  levels  using  T  (700)  as  the  forcing-level 
temperature. 

T  =  -268.4072  +  0.1473T3  +  1.0  X± 

X,  =  0.555Z.,  +  0.532Z„  -  0.521Z_  +  0.061ZC  -  0.035Znc  +  0.081Z,, 

1  1         I  3  5         15         16 

-0.059Zlg 

T2  =  -285.3345  -  0.0541T3  +  1.0  X2 

X„  =  -Q.014Z,  -  0.016Z„  +  0.319Z-  +  Q.324Z.  -  0.020Z.,  +  0.008Zn 

2  12         3         4         7         9 

T.  =  -293.3438  -  0.0755To  +  1.0  X. 
4  3        4 

X.  =  0.106Zri  +  0.109Zin  +  0.101Z,,  +  0.095Zin  -  0.036Z1o 

4  9         10         11         12         13 

T  =  -277.0920  +  0.0092T  +  1.0  X 

Xc  =  -0.041Z1o  +  0.171Z,.  +  0.094Z1t  -  0.014Z-, 

5  13         14         15         16 

T,  =  -244.0715  +  0.0905To  +  1.0  X, 

6  3        6 

X,  =  -0.026Z.  -  0.017Z.,  -  0.015Znc  +  0.059Z,,  +  0.201Z...,  -  0.026Zno 
6  4  7  15  16  17  18 

T?  =  -262.7598  +  0.0306T3  +  1.0  Xy 

X,  =  -0.020Z.  -  0.037Zn^  +  0.145Zlo  +  0.114Zin  -  0.018Z__ 

Tg  =  -272.9868  +  0.0065T3  +  1.0  Xg 

Xg  =  -0.003Z12  -  0.026Z   +  0.099Z2Q  +  0.105Z   -  0.047Z22 

Tg  =  -254.3952  +  0.0447T3  +  1.0  Xg 

Xn  =  -0.022Z.  -  0.008Z-  -  0.005Z, _  -  0.029Zo1  +  0.308Zo„ 
9  4         9         15  21  22 

X-.  is  specified  in  terms  of  heights  z  as 

X1  =  0.555(zg50-z1000)  +  0.532(zgoo-z950)  -  0. 521(zg50-zgoo) 

+0.061(z750-z800)  -  0.035(z250-z300)  +  O.OoKz^-z^) 

-0.059(z150-z175). 
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TABLE  5.   Full-sample  regression  equations,  with  associated  thickness- 
variables  at  Resolute,  NWT,  for  the  specification  of  the  January- 
temperatures  at  indicated  levels  using  T  C5Q0)  as  the  forcing-level 
temperature. 

Tn  =  -318.9878  -  0.0034T.  +  1.0  X, 
1  4        1 

X,    =     0.410Z.,   +  0.426Zo   -   0.392Zo  +  0.233Z.    +  0.010Z1C  +  0.068Z,o 

1  1  2  3  4  15  18 

-0.056Z22 

T0  =  -269.0266  +  0.0200T.    +  1.0  X- 

2  4  2 

X_  =   -0.022Z,    -  0.009Z„  +  0.325Z-  +  0.284Z.    -  0.008Z, 

2  12  3  4  6 

T0  =  -289.0706  =   0.0791T.   +  1.0  X. 

3  4  3 

X„   =  -0.090Z.   +  0.169Zq  +  0.133Z.,  +  0.124Zo  +  0.170Z,   -  0.011Z- 
3  4  5/869 

Tc   =   -262.3394  +  0.0421T.   +  1.0  Xc 
5  4  5 

Xc  =  -0.007Z,    -  0.011Z1o   -   0.034Z1o  +  0.159Z,.   +  0.10QZ1C   -  0.017Z,, 

5  1  12  13  14  15  16 

T,   =   -272.8584  +  0.0232T.    +  1.0  X, 

6  4  6 

X,   =   -0.011Z1C   +  0.064Z.,   +  0.206Z..,   -   0.035Zlo 

6  15  16  17  18 

T-.   =   -264.4431  +  0.0184T.   +  1.0  X_ 

7  4  7 

X.,  =   -0.01lZ-,„    -   0.032Z.,-.  +  0.133Zno  +  0.122Zin    -  0.018Zon 

/        iz       i/       io       iy       zu 

T0  =  -270.4817  +  0.0114T.  +  1.0  X0 
6  4         b 

X0  =  -0.005Z..  +  0.020Znn  +  0.083Z„n  +  0.121Zon  -  0.058Zoo 
o        /        iy        zu        zi        az 

T„  =  -269.8740  -  0.0072T.  +  1.0  Xn 
9  4       9 

Xn  =  -0.005Z„  +  0.004Z1/  -  0.005Z-,,  -  0.040Zon  +  0.325Zo„ 
9  z         14         16         21         11 
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TABLE  6.   Full-sample  regression  equations,  with  associated  thickness- 
variables  at  Thule,  Greenland,  for  the  specification  of  the  January 
temperatures  at  indicated  levels  using  T  C3Q0)  as  the  forcing-level 
temperature. 

T  =  -257.2419  -  0.0353T5  +  1.0  X1 

Xn    =     0.534Z-,   +  0.554Z„   -  0.364Z-   -  0.037Zo   -  0.041Ziri  -  0.022Zn, 

1  1         2         J         9         10         16 

+0.043Zlg 
T  =  -272.4780  +  0.0045T5  +  1.0  X£ 

X0  =  -0.019Zo  +  0.338Zo  +  Q.315Z.  -  0.020ZC  -  0.03QZ, 

2  2         3         4         5         6 

T  =  -272.8140  +  0.0Q93T5  +  l.Q  X3 

X„  =  0.057Zo  -  0.118Z.  +  0.166Z.  +  0.162Z,  +  0.108Z.,  +  0.102Zo 

3  3  4  5  6/8 

T.   =   -257.5325  +  0.0588TC  +  1.0  X. 

4  5  4 

X.    =     0.074Z„  +  0.073Znn  +  0.119Zn1    +  0.125Zlo   -   O.Q42Z1„   -  0.019Z-. 
4  9  10  11  12  13  14 

T,   =   -246.8040  +  0.0773T.  +   1.0  X, 
6  5  6 

X,   =   -0.008Z,    -   0.035Zn[.  +  0.059Z..,   +  0.198Z...,   -   0.013Z1o 
6  1  15  16  1/  18 

T     =  -270.2495  +  0.0083T5  +  1.0   X? 

X-,  =  -0.009Z,    +  0.006Z^   -   0.038Z^  +  0.156Z1o  +  0.081Z1o 
/  i  /  1/  io  xy 

Tg  =   -275.9387  +  0.0118T5  +  1.0   Xg 

X0  =  -0.006Z„  +  0.012Z.  +  0.061Zon  +  0.114Zon  -  0.036Z„o 

8  2  4  20  21  11 

Tg  =  -283.4604  -  0.0250T5  +  1.0  Xg 

X„  =  -0.006Z,  +  0.011Ziri  +  0.012Zn,  -  0.027Zn„  +  0.277Z„o 

9  1         10         16         19         22 
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TABLE  7.   Full-sample  regression  equations,  with  associated  thickness- 
variables  at  Thule,  Greenland,  for  the  specification  of  the  January 
temperatures  at  indicated  levels  using  T  (200)  as  the  forcing-level 
temperature. 

T.,  =  -331.3064  -  0.2659T,  +  1.0  X.. 

1  6        1 

X  =  0.526Z1  +  0.543Z2  -  0.346Z3  -  0.035Zg  -  0.039Z1Q  +  Q.Q86Zlg 
-0.156Z22 

T0  =  -272.2056  -  0.0038T,  +  1.0  X„ 

2  fa         z 

X0  =  -0.019Zo  +  Q.337Z-  +  Q.314Z,  -  Q.Q20ZC  -  0.Q29Z, 

2  2         3         4         5         6 

T0  =  -273.0007  +  0.Q049T,  +  1.0  X„ 

3  6         3 

X0  =  0.057Z„  -  0.118Z.  +  0.165ZC  +  0.162Z,  +  0.109Z^  +  0.102Zo 

3  3  4  5  6  7  8 

T,   =   -269.6765  +  0.0154T,   +   1.0  X. 

4  6  4 

X.    =     0.074Z„  +  0.074Z,„  +  0.118Zni    +  0.124Z,„   -   0.045Zno   -  0.006Zn. 

4  9         10         11         12         13         14 

Tc  =  -264.6812  +  0.0508T,  +  1.0  Xc 

5  6        5 

Xc  =  0.009Zr  -  0.055Z1o  +  0.179Z1/  +  0.091Z,.  -  0.019Z., 
5  6         13         14         15         16 

T-,  =  -195.1660  +  0.2897T,  +  1.0  X., 
7  6        7 

X^  =      0.010Z,.   -   0.021Z,,    -  0.091Z,.,  +  0.161Z1Q  +  0.081Z1o 
7  15  16  1/  lo  19 

T0  =   -283.8765   -   0.0262T.   +  1.0  X0 
o  bo 

Xg  =      0.00GZ     +  0.008Z       +  0.060Z2Q  +  0.114Z        -  0.036Z22 

Tn   =   -280.3738   -   0.0026T,   +   1.0   Xn 
y  oy 

X  =   0.009Zg  +  0.010Z16  -  0.026Z   +  0.277Z20 
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TABLE  8.   Full-sample  regression  equations,  with,  associated  thickness- 
variables  at  Alert,  NWT,  for  the  specification  of  the  July  tempera- 
tures at  indicated  levels  using  T  (200)  as  the  forcing-level 
temperature. 

T,  =  -273.6511  -  0.3203T,  +  1.0  Xn 
1  6         1 

X,  =  0.425Z  +  0.373Z2  -  0.240Z3  -  0.033Zg  -  0.027Z  +  0.029Z12 
+0.032Z16 


T  =  -304.2915  -  0.0173T6  +  1.0  X2 


X0  =  0.026Z,  +  0.024Z„  +  0.142Z„  +  Q.399Z.  -  0.009Z^  -  0.007Znc 

2  1         2         3         4         7         15 

+0.036Z1o 
lo 

T0  =  -273.5454  -  0.0034T,  +  1.0  X, 

3  6        3 

Xn  =  0.034Zo  -  0.132Z.  +  0.164ZC  +  0.172Z,  +  0.127Z^  +  0.114Zo 

3  3  4  5  6  7  8 

-0.007Z 

T.    =   -278.9707  +  0.0340T,   +  1.0  X. 

4  6  4 

X.    =     0.006Z„  +  0.078Zn  +  0.082Zin  +  0.120Znn   +  0.112Zn_   -  0.053Zno 

4  Z  y  11)  11  12  13 

+0.010Z-,  , 
14 

Tc  =  -248.0954  +  0.1898T,  +  1.0  Xc 

5  6        5 

Xc  =  -0.011ZC  +  0.047Z,,  +  0.032Z10  -  0.174Z, „  +  0.264Zn.  +  0.082Z1C 
5  5         11         12         13         14         15 

-0.040Z, , 
16 

T,  =  -258.6167  +  0.0495T,  +  1.0  X., 
/  6        7 

X..  =  0.015Z.  -  0.010Zin  -  0.004Z, _  -  0.047Zn^  +  0.155Z1o  +  0.101Z10 
/  4  11  15  1  /  lo  iy 

-  0.014Z2Q 
T0  =  -272.3201  +  0.0112T,  +  1.0  X0 

O  DO 

X0  =  -0.006Z_  +  0.010Z.  -  0.004Z. .  -  0.004Z.,  +  0.061Z„n  +  0.134Zol 

8  2         4  14  16  20  21 

-0.063Z22 

Tn  =  -273.5842  +  0.0138T,  +  1.0  Xn 

9  6        9 

X„  =  -0.006Z..  -  0.010Z„  +  0.022Z_  -  0.005Ziri  +  0.008Z._  -  0.Q13Z-- 
y        j       /       3       iu       13       zi 

+0.271Z22 
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TABLE  9.   Full-sample  regression  equation,  with  associated  thickness- 
variables  at  Alert,  NWT,  for  the  specification  of  the  July  tempera- 
tures at  indicated  levels  using  T  (150)  as  the  forcing-level 
temperature. 

T  =  -207.5527  -  0.0898T?  +  1.0  X1 

X,  =  0.424Zn  +  0.374Zo  -  0.248Zo  +  0.024Z.  -  0.045Z.  -  0.035Zn 
i        l  i  5  o       o       y 

+0.022Z12 
T2  =  -329.3796  -  0.1059T?  +  1.0  X£ 

X„  =  0.026Z,  +  0.024Z.  +  0.149Z.  +  0.387Z.  -  0.008Zn.  +  0.055Zn^ 

2  1         2         3         4         14         17 

T  =  -274.1506  -  0.0104T  +  1.0  X 

X0  =  0.033Zo  -  0.127Z.  +  0.160ZC  +  0.168Z£  +  0.125Z.,  +  0.112Zo 

3  3  4  5  6  7  8 

T.    =  -271.5181  +  0.03107V    +   1.0  X. 

4  7  4 

X.    =     0.080Zo  +  0.084Ziri  +  0.112Zn1   +  0.104Z--  -   0.037Z1o 
h  y  xu  -Li.  j-Z  li 

T5   =  -276.9827  +  0.1162T     +  1.0  X 

X,  =  0.047Zni  +  0.032Z1o  -  0.175Zlo  +  0.264Zn.  +  0.081Znc  -  0.025Z1£ 

5  11         12         13         14         15         16 

T,  =  -259.7761  +  0. 1023T.,  +  1.0  X, 

6  7        6 

X,  =  0.007Z_  +  0.038Zn,  +  0.184Z,,  -0.023Znn  +  0.016Zon 
6  1         16         17        19         21 

T0  =  -285.4023  -  0.0316T-  +  1.0  X0 

X0  =  -0.005Zo  +   0.010Z.    -   0.002Z1C  +  0.062Z^  +   0.136Zon    -   0.062Zoo 
o  Z  4  15  20  ZL  22 

Tg   =  -273.5064  +  0.0135T     +   1.0  X 

X_   =  -0.007Z_  +  0.013Z.   -   0.004Znn  +  0.007Z1o   -   0.012Zol   +   0.270Zoo 
y  Z  J  11)  13  2L  zz 
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TABLE  10.   Full-sample  regression  equation,  with,  associated  thickness- 
variables  at  Thule,  Greenland ,  for  the  specification  of  the  July 
temperatures  at  indicated  levels  using  T  (1Q0)  as  the  forcing-level 
temperature. 

T.  =  -215.4731  +  0.3238TQ  +  1.0  X1 
1  o        1 

X,  =  0.360Z-.  +  0.374Z„  -  0.075Zo  -  0.034Z.  -  0.041Z,  +  0.020Znc 

1  1         l  3         4         6         15 

-0.021Z,, 
16 

T2  =  -303.4197  -  0.0272Tg  +  l.Q  Xg 

X„  =  0.047Z,  +  0.048Z„  +  0.105Zo  +  0.401Z.  -  0.012ZC  -  0.004Zo 

2  12  3         4         5  8 

+0.018Z 

T,  =  -319.8106  -  0.137QTo  +  1.0  X„ 
j  o         J 

X„  =  -0.081Z.  +  0.164ZC  +  0.156Z,  +  O.lllZ^  +  0.116Zo  +  0.047Zon 

3  4         5         6         7         8         21 

-0.031Z22 

T,  =  -276.6890  +  0.0313To  +  1.0  X. 

4  8        4 

X.  =  0.011Zo  +  0.092Zn  +  0.094Zin  +  0.096Z__  +  0.096Z, „  -  0.033Z1o 

4  2  9  10         11  12  13 

Tc  =  -266.3704  +  0.0764To  +  1.0  Xc 

5  o         d 

Xc  =  0.014Z.,  -  0.073Zno  +  0.214Z-,  +  0.070Z,C  -  0.016Z,, 

5  7         13         14         15         16 

TV  =  -255.2715  +  0.1766To  +  1.0  X, 

6  8        6 

X,  =  0.019Zn.  -  0.035Znc:  +  0.061Z..,  +  0.123Z,.,  +  0.103Zno  -  0.044Z„„ 

6  14         15         16         17         18         20 

T.,  =  -184.3870  +  0.0909To  +  1.0  X^ 

X.,  =  -0.030Z,  -  0.050Z-,  +  0.031Z...  -  0.014Z.,,  +  0.068Z^  +  0.008Z1o 

7  6         9         11         16         17         18 

+0.021Z 
Tn  =  -257.1123  -  0.020TQ  +  1.0  Xn 

y  o      y 

Xn   =     0.011Z.   +  0.002Znn    -   0.008Zlo   -   0.003Z...    +  0.24lZoo 
y  4  11  13  14  11 
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B.   RELATED  STATISTICAL  PARAMETERS 

Available  as  outputs  from  the  BIMED  02R  program  are  several  related 
statistical  parameters  such  as: 

a.  multiple  R  after  r  steps 

b.  standard  error  of  estimate,  S.E. 

c.  mean  value,  T 

d.  standard  deviation,  a 

e.  F-statistic  after  r  steps,  F^ 

f.  percent  reduced  variance,  Rz 


r 

I2. 
r 


2 
The  significance  of  R  in  relation  to  S.E.  is  given  by  (Crow  et  al. , 


1955)  as 


where 


(S.E.)2  =  a2((N-l)/(N-r-l))(l-R2) 


2     N     -  2 
o   =   E  (T.-Tr  /  N.  (3) 

i=l  X 
It  should  be  noted  that  S.E.  is  a  measure  of  the  error  in  specification 
of  T   after  application  of  the  regression  equation  estimator  of 

form  (2). 

2 
In  (3) ,  a   is  the  variance  of  the  temperature  sample  and 

N   is  the  sample  size 

i   is  the  sample- element  index 

r   is  the  number  of  predictors  selected 

R  is  the  multiple  correlation  coefficient  after  r  steps.   It 

2 
follows  that  the  percent  unexplained  variance  (1  -  R  ) ,  for  sample  sizes 

of  213  to  244,  as  contained  here,  is  closely  approximated  by 

(S.E.  /  a)2  =  1  -  R2. 
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2 
The  percent  explained  variance  R  ,  also  called  the  percent  reduced 

variance,  is  then  given  by 

R2  =  1  -  (S.E.  /  a)2.  (4) 

2 
The  mean  value  of  R  ,  found  in  the  prediction  of  the  1000  mb  tempera- 
tures using  the  full-data  set,  was  0.8478  and  0.7480  with  the  10% 

2 
extreme-data  set.   For  all  levels  above  1000  mb ,  the  pooled  mean  R  for 

the  full-data  set  was  0.9672  and  0.9596  for  the  10%  extreme-data  set. 

The  scores  at  1000  mb  are  somewhat  lower  than  those  of  the  upper  levels 

and  this  is  to  be  expected,  in  that  the  surface  variabilities  bring 

periodic  and  random-effect  terms  into  the  specification  of  the 

temperature  at  1000  mb. 
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IV.   ANALYTICAL  PROCEDURE  USING  THE  REGRESSION  TECHNIQUE 

A.   RESULTS  OF  THE  REGRESSIONS 

Part  (a)  of  Tables  11  through  17  lists  the  mean,  standard  deviation, 
multiple  correlation  coefficient  and  standard  error  of  estimate  for  the 
full-data  sample  of  each  station  and  at  all  mandatory  levels  except  the 
forcing  level.   These  statistics  were  computed  after  application  of  the 
appropriate  regression  equations  listed  in  Tables  4  through  10.   Part 
(b)  of  Tables  11  through  17  gives  similar  results,  after  application  of 
the  identical  regression  equations  used  in  part  (a) ,  to  a  regression 
analysis  of  the  data  set  which  contains  only  the  "10%  extreme"  set  of 
forcing -level  sample  atmospheres.   This  means  that  T(P  )  was  chosen  if 
it  lay  within  the  particular  10%  Gaussian  distribution  extreme  under 
consideration.   The  remainder  of  the  mandatory  levels  of  the  extreme 
sounding,  T (. 10) ,  were  then  subjected  to  the  same  regression  analysis 
as  the  full-data  samples  of  part  (a) .   This  nominal  set  of  extreme 
atmospheres  was  selected  by  choosing  those  soundings  for  which 
T(P,)  >  T  (.10)  for  the  warm  extremes,  and  T  (P  )  <_  T  (.10)  for  the  cold 
extremes.   These  extreme  10%  values  T  (.10)  are  defined  by  the  Gaussian 
distribution  limits  listed  in  Eq.  (1). 

Part  (c)  lists  the  mean  temperature  profile  defined  at  mandatory 
levels  which  correspond  to  the  set  of  soundings  having  T (P  )  at  the  1% 
extreme.   There  were  two  ways  of  selecting  the  1%  extreme  cases:   (i) 
application  of  an  a  priori  test  to  determine  the  values  T(P  )  which 
fell  within  the  required  range  indicated  by  the  Gaussian  distribution; 
(ii)  if  there  were  no  values  of  T (P  )  fitting  this  constraint,  then 
either  the  two  or  three  coldest  (for  cold  extremes)  or  warmest  values, 
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were  chosen  a_  posteriori  from  the  regression-determined  1%  set  as  repre- 
senting the  set  of  soundings  corresponding  to  the  1%  extreme  T(P  )- 
values.   The  corresponding  set  of  T(P  )  values  were  then  simply  averaged 
and  considered  to  represent  the  expected  vertical  profile  having  T(P  ) 
in  the  1%  extreme. 


B.   THE  TEST  FOR  SIGNIFICANT  SOUNDING  DIFFERENCES 

In  view  of  the  small  sample  sizes  in  parts  (b)  and  (c)  of 
Tables  11,...,  17,  it  was  not  feasible  to  test  for  significant  differ- 
ence by  random  sampling  of  the  10%  set  relative  to  the  1%  set  of  T(P  ) 
generated  soundings.   It  was  decided  therefore  to  test  for  significant 
differences  in  defining  the  T  values  predicted  in  each  part  (b)  of 
Tables  11,...,  17  against  those  T  values  of  part  (a)  obtained  by  the 
identical  prediction  process,  but  using  the  full-data  sample:   that  is, 
in  case  (a)  there  were  no  a_  priori  limiting  conditions  set  up  regarding 
a  particular  value  of  T(P  ). 

In  order  to  check  the  latter  hypothesis,  a  Student's  t-test  was 
performed  to  test  whether  the  10%  samples  of  prediction  means  of  part  (b) 
were  significantly  different  from  those  of  the  full-data  sample  means 
of  part  (a).   The  pooled  t-statistic  relating  to  the  T  prediction 
differences  of  the  part  (a)  and  part  (b)  means  may  be  formulated  after 
Brooks  and  Carruthers  (1953)  with  (N+n-(r+2))  degrees  of  freedom,  after 
computing  the  two  sets  of  predictive  mean  values  and  their  standard 
errors.   The  t-statistic  is  formulated  as  follows: 


t(N+n-(r+2))  =  (T^-T^) 


(N-9)o*+  (n-2)o*    X*/l+l  ^ 
(N+n-r-2)      7  \N   " 


(5) 
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Here  the  T's  are  the  part  (a)  and  the  part  (b)  sample  temperature  means 

T   (at  the  mandatory  levels)  based  upon  identical  prediction  equations. 

o     is  the  part  (a)  standard  error  of  estimate,  and  a       in  is  the  part 
E  b , • 1U 

(b)  standard  error  of  estimate.   The  number  of  degrees  of  freedom  asso- 
ciated with  the  two  a^,  values  are  (N-r-1)  and  (n-2)  ,  respectively,  in 
parts  (a)  and  (b)  with 

N  =  the  full-sample  size 

n  =  the  nominal  10%  extreme-sample  size 

r  =  the  number  of  predictors,  r  =  8. 
The  pooled  t-statistic  may  be  computed  from  the  horizontally  adjacent 
sections  of  Tables  11,...,  17  and  with  the  proper  number  of  (N+n- (r+2)) 
degrees  of  freedom  for  t.   Tables  18  and  19  of  Appendix  1  list  the 
resulting  values  at  mandatory  levels  for  each  station  under  consideration. 

These  t-values  were  tested  against  critical  t   ..-values  (e.g., 
Crow  et  al,  1960)  with  the  same  number  of  degrees  of  freedom  used  in 
(5).   This  number  of  degrees  of  freedom  usually  fell  into  the  tabular 
class  for  t-values  corresponding  to  an  infinite  number  of  degrees  of 
freedom.   Thus  for  example  t    ,  with  an  infinite  number  of  degrees  of 
freedom,  is  t  =  2.326.   Tables  18  and  19  show  one  value  below  the  99% 
confidence  limit  specified  by  t    =  2.326.   That  value,  corresponding 
to  the  1000  mb  level  of  Thule  (with  P   =  300  mb) ,  still  yields  a  con- 
fidence limit  of  92%.   The  random  chance  of  the  T   specifications  of 

M 

parts  (a)  and  (b)  of  Tables  11,...,  17  being  drawn  from  the  same  popu- 

12 
lation  if  t  =  7.0  is  approximately  one  chance  in  5  X  10   and  is  much 

smaller  for  t  >  7.0.   Most  of  the  t-values  in  Tables  18  and  19 

exceeded  t  =  7.0  with  exceptions  principally  at  P  =  1000  rab. 
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C.   THE  TYPE-(C)  PROFILE  CLASSIFICATION 

As  noted  previously,  it  was  not  feasible  to  employ  a  t-test  utiliz- 
ing the  extremely  small  sample  size  of  part  (c)  of  Tables  11,...,  17. 
However,  it  is  useful  to  explain  physically  the  pattern  of  trends  in 
the  sample  means  T   of  part  (c)  relative  to  those  of  (b) .   The  code 
symbol  "D"  for  "direct",  used  in  part  (c) ,  means  that  T  (c)  is  dis- 
placed directly  in  the  same  sense  from  T  (b)  as  T  (b)  is  from  T  (a)  at 
the  level  under  consideration.   The  code  symbol  "M"  for  "mixed"  means 
that  |t(c)  -  T  (a) |  <   lTM(b)  ~  TM(a) I  so  that  a  consistent  displace- 
ment in  the  temperature  extreme  does  not  occur  in  progressing  across  the 

columns  of  T  -values.   This  code  symbolism  was  used  for  all  cases  con- 
M 

sidered  and  is  summarized  under  part  (c)  of  Tables  11,...,  17.   To 
some  extent,  the  nature  of  these  displacements,  whether  direct  or  mixed, 
is  indicated  also  by  Fig.  2,...,  Fig.  8,  which  show  only  the  comparison 
of  the  T  (c)  and  T  (a)  profiles. 


38 


V.   RESULTS  SUMMARY  FOR  INDIVIDUAL  STATIONS 

A.   WINTER  CASES 

This  section  will  consider  those  cases  investigated  as  cold  ex- 
tremes and  summarized  in  Tables  11,...,  14  and  Figs.  2,...,  5. 

1.   Hall  Beach,  NWT;  T(Pj)  =  ^(700) 

Table  11(a)  shows  the  results  of  the  stepwise  regression 
applied  to  the  specification  of  the  January  full-sample  temperatures 
(T   =  T  ,  T  ,  T  ,  T  ,  Tfi,  T  ,  T   and  T  )  at  the  mandatory  levels  of 
Hall  Beach  by  an  equation  of  form  (2) .   An  eight-predictor  equation  for 
T  ,  at  each  level,  was  generated  by  use  of  the  BIMED  02R  program  with 
the  results  shown.   Very  high  multiple  correlation  coefficients  were 
obtained,  indicating  a  functional  relationship  for  T   at  all  levels. 
Table  11(b)  lists  analagous  results  for  the  data  sample  corresponding 
to  the  10%  cold  extreme  at  the  forcing  level,  with  T  (b)  being  deter- 
mined by  the  same  specification  equation  used  in  part  (a) .   Part  (c) 
presents  the  resultant  estimated  mandatory-level  means  T  (c)  derived 
from  the  five  cases  which  fell  into  the  1%  extreme  forcing  level  sub- 
set of  part  (b) .   Part  (c)  then,  gives  the  expected  vertical  profile 
corresponding  to  the  1%  extreme  at  T(P  )  =  T(700).   It  is  suggested 
that  the  mean  T  (c)  profile,  relative  to  that  of  T  (b)  ,  should  be  phys- 
ically as  well  as  statistically  explainable  from  the  conditions  imposed 
by  the  choices  of  P  .   At  Hall  Beach,  the  latter  selection  at  700  mb 
was  consistent  with  the  existence  of  an  extreme-cold  troposphere, 
T  (c)  <  T  (b) ,  and  an  extreme-warm  stratosphere  for  winter.   The  only 
case  of  a  mixed  classification  M,  occurs  for  T  (500) ,  where  according 
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FIG.  2.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Hall  Beach,  NOT,  corresponding 
to  the  set  of  1%  cold  extreme  occurrences  of  T  (700)  .   The  thin 
solid  line  depicts  the  January  mean  Hall  Beach  vertical  sounding 
(1967-70). 
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to  Fig.  2,  the  crossover  between  the  mean  and  extreme  soundings  occurs. 
The  nature  of  the  TM(c)  profile  is  then  consistent  with  the  displacement 
of  the  TM(b)  -  TM(a)  profiles  except  at  the  crossover  point  near  the 
tropopause. 

2.   Resolute,  NWT;  T (P  )  =  T  (500) 

In  the  case  of  Resolute,  Table  12(c)  shows  a  somewhat  different 
result  from  that  discussed  above  for  Hall  Beach.   At  Resolute  we  now 
start  with  the  premise  of  an  initially  specified  very  cold 
T  (c)  <  T    (500).   This  level  is  excluded  from  being  an  effective 
crossover  point  in  the  results  of  Fig.  3  except  in  the  region  of  P=300  mb 
where  a  slight  tendency  toward  a  warm-stratosphere  does  occur.   Compar- 
ing the  two  cases,  the  selection  of  an  extreme  cold  T  (500)  dictates  the 
presence  of  a  cold  stratosphere  coupled  with  a  cold  troposphere  below. 

For  nearly  all  levels  at  Resolute,  the  T  (c)  -  T  (a)  trend  bears  out 

M        M 

that  indicated  by  T  (b)  -  T  (a)  in  the  sense  of  maintaining  the  T  (c) 
atmosphere  relatively  cold  at  virtually  all  levels.  In  contrast  with 
the  situation  at  Hall  Beach,  the  Resolute  type  (c)  atmosphere  more 
nearly  resembles  the  case  of  a  cold  tropospheric  vortex  with  an  over- 
lying stratospheric  cold  vortex,  whereas  the  Hall  Beach  type  (c)  case 
seems  to  occur  with  the  warm  stratospheric  anticyclone  aloft. 


3.   Thule,  Greenland;  T (P  )  =  T  (300) 

This  case  is  graphically  depicted  in  Fig.  4,  which  shows  that 
the  _a  priori  condition  of  a  1%  cold  temperature  extreme  at  P  =300  mb  at 
Thule  is  associated  with  a  cold  stratosphere  as  well  as  a  cold  tropo- 
sphere at  all  levels  above  and  below  P  =300  mb .   The  extreme  vertical 
profile  for  this  regime  0(1'  )  =  T(300))  at  Thule  resembles  that  at 
Resolute,  where  there  was  also  evidence  of  a  similar  coupling  between 
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FIG.  3.  The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Resolute,  NWT,  corresponding  to 
the  set  of  1%  cold  extreme  occurrences  of  TT(500).  The  thin  solid 
line  depicts  the  January  mean  Resolute  sounding  (1967-70). 
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FIC.  4.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Thule ,  Greenland,  corresponding 
to  the  1%  cold  extreme  occurrences  of  T  (300) .   The  thin  solid  line 
depicts  the  January  mean  Thule  sounding  (1967-70) . 
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FIG.  5.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Thule,  Greenland,  corresponding 
to  the  1%  cold  extreme  occurrences  of  T  (200) .   The  thin  solid  line 
depicts  the  January  mean  Thule  sounding  (1967-70). 
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cold  stratospheric  temperatures  in  case  (c)  and  a  cold-extreme  upper 
troposphere.   Now,  however,  Fig.  4  shows  no  crossover  point  between  the 
mean  and  extreme  profiles.   Likewise  Table  13  shows  the. direct  transi- 
tion D-type  classification  applicable  at  all  levels  of  the  mean 
temperature  profile  T  (c) . 

4.  Thule,  Greenland;  T(P  )  =  T  (200) 

This  case  differs  from  the  previous  case  (T(P  )  =  T  (300)),  in 
that  the  latter  was  associated  by  regression  development  with  T  (200)  = 
-62. 83°C,'  whereas  the  coldest  1%  input  temperature  T  (200)  was  -74.70°C. 
The  imposition  of  the  latter  condition  into  the  regression  analysis  at 
Thule  (T  (200))  is  associated  with  the  coldest  stratospheric  regime 
encountered  during  the  entire  study  of  this  paper.   With  T  (200)  = 
-74.70  C,  all  regression-generated  stratospheric  temperatures  are  typi- 
cal of  those  in  the  cold  Arctic  vortex  at  and  above  200  mb .   However 
Fig.  5,  shows  that  in  this  case  there  is  a  crossover  in  the  T  (a) ,  TM(°) 
profiles.   Furthermore  by  comparison  of  Table  14(c)  with  Tables  11(c), 
12(c)  and  13(c),  it  is  clear  that  this  case  corresponds  to  the  warmest 
of  the  four  winter  tropospheres  discussed  here. 

5.  Combined  Winter  Regimes,  Summary 

The  analysis  of  the  four  sets  of  type  (c)  winter  extremes  has 
indicated  that  the  300  mb  level  is  a  key  indicator  in  the  type  of  atmo- 
sphere to  be  expected  in  the  Canadian  Arctic  region.   In  summary, 
these  regimes  have  been  identified  by  the  four  stratifications: 

(a)  Hall  Beach  T  (700) ,  cold  extreme 

(b)  Resolute   T  (500) ,  cold  extreme 

(c)  Thule     T .(300) ,  cold  extreme 

(d)  Thule      T  (200),  cold  extreme 
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Case  (a)  corresponds  to  an  abnormally  cold  troposphere  with  a 
resultant  abnormally  warm  stratosphere  (Fig.  2).   Cases  (b)  and  (c) 
correspond  jointly  to  cold  stratospheres  and  cold  tropospheres  (Figs.  3 
and  A),  but  neither  cold  layer  exists  at  a  particular  extreme.   Finally 
case  (d)  is  characterized  by  a  record  cold  stratosphere  and  a  near- 
record  warm  lower  troposphere.   Cases  (a)  and  (d)  show  negative  corre- 
lations of  temperature  across  the  tropopause,  a  property  documented  at 
mid-latitude  stations  by  Cole  and  Nee  (1965) .   However  cases  (b)  and 
(c)  have  a  consistent  positive  correlation  throughout  the  atmosphere 
and  have  yet  to  be  documented  fully. 

B.   SUMMER  CASES 

This  section  will  consider  those  cases  as  having  warm  extremes  at 
forcing  levels  and  are  summarized  in  Tables  15,  16  and  17  and  Figs.  6, 
7  and  8.   The  same  stepwise  regression  procedure  was  applied  to  the 
July  data  of  the  warm  extreme  stations  as  was  applied  to  the  January 
data  of  the  cold  extreme  stations.   The  results  are  presented  below. 

1.   Alert,  NWT;  T(PJ)  =  T  (200) 

Table  15  lists  the  statistical  results  of  the  analysis  of  Alert 
with  T(P  )  =  T  (200.   Consistently  high  multiple  correlations  are 
shown  indicating  a  high  degree  of  confidence  in  the  predicted  T  values 
including  those  of  part  (c) .   This  confidence  is  reinforced  by  the  con- 
sistent direct  transition  D-type  classification  applicable  at  most 
levels  of  the  mean  temperature  profile  T.,(c).   Mixed  tendencies  occur 
at  100  mb  and  70  mb  but  are  characterized  by  extremely  small  tempera- 
ture differentials  (less  than  0.5  C)  in  the  T  (c)  -  T  (h)  trends.   A 

crossover  of  the  l'(c)  extreme  and  T.,(a)  profiles  exists  between  500  mb 
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FIG.  6.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Alert,  NWT,  corresponding  to  the 
1%  warm  extreme  occurrence  of  T  (200)  .   The  thin  solid  line  depicts 
the  July  mean  Alert  sounding  (1967-70) . 
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and  300  mb  (Fig.  6),  indicating  an  extreme  warm  stratosphere  overlying 
a  relatively  cool  troposphere.   This  phenomenon  is  also  characteristic 
of  mid-latitude  negative  inter-level  temperature  correlations  shown  by 
Cole  and  Nee  (1965)  to  exist  near  the  level  of  the  tropopause  during 
periods  of  cyclone-anticyclone  variability. 

2.  Alert,  NWT;  T(P  )  =  T  (150) 

Table  16  corresponds  to  the  analysis  of  Alert  with  T(P  )  = 
T  (150)  .   Very  high  multiple  correlation  coefficients  in  both  parts  (a) 
and  (b)  are  again  shown  yielding  a  resultant  high  degree  of  confidence 
in  the  diagnostic  procedure  of  isolating  "extreme"  atmospheres. 
Figure  7  indicates  that  a  marked  similarity  exists  between  the  T  (a) 
and  T  (c)  profile  pairs  and  the  corresponding  pair  for  Alert  of  T  (200) . 
Table  16  indicates  that  the  T  (c)  regression-generated  atmosphere  is 
displaced  primarily  in  a  "direct"  sense  relative  to  the  T  (b)  atmo- 
sphere.  It  should  also  be  noted  that  the  forcing  extreme  T  (150)  at 
the  nominal  10%  and  1%  values  differ  only  by  about  1  C.   The  comparison 
between  T  (c)  and  T  (b)  is  generally  of  the  same  character  as  that  of 
TM(c)  relative  to  T  (a)  ,  (Fig.  7). 

The  same  conditions  of  warm  stratosphere  and  relatively  cool 
troposphere,  as  well  as  the  crossover  between  500  mb  and  300  mb,  sug- 
gest a  similar  comparison  with  regard  to  this  climatic  regime  as  with 
the  preceding  one  (Alert  T  (200) ) . 

3.  Thule,  Greenland;  T(P})  =  TjQOO) 

The  final  table  of  regression  statistics,  Table  17,  corresponds 
to  the  analysis  of  Thule  at  T(P  )  =  T  (100)  .   Once  again  there  are  very 
high  multiple  correlation  coefficients  in  both  parts  (a)  and  (l>).   In 
addition  the  atmospheric  profiles,  both  in  the  mean  and  1%  extn 
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FIG.  7.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Alert,  NWT ,  corresponding  to  the 
1%  warm  extreme  occurrence  of  T  (150) .   The  thin  solid  line  depicts 
the  July  mean  Alert  sounding  (1967-70) . 
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FIG.  8.  The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Thule,  Greenland,  corresponding 
to  the  1%  warm  extreme  occurrence  of  T  (100) .  The  thin  solid  line 
depicts  the  July  mean  Thule  sounding  (1967-70) . 
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(Fig.  8),  were  nearly  identical  to  the  corresponding  profile  pairs  in 
the  previous  two  warm  cases  (Figs.  6  and  7).   This  indicates  a  strong 
tendency  of  the  synoptic  processes  acting  at  Alert  in  July  to  also  be 
acting  at  Thule  during  the  same  period. 
4.   Combined  Summer  Regimes,  Summary 

Each  of  the  July  cases  is  based  upon  the  assumption  of  a  warm 
extreme  forcing- level  temperature  existing  at  or  above  200  mb.   In  each 
case,  there  proved  to  be  a  concurrent  warm  temperature  at  300  mb,  that 
is,  stratospheric  warming  above  200  mb  has  extended  downward  to  300  mb 
in  July.   At  tropospheric  levels,  however,  the  mean  temperatures  of 
parts,  (b)  and  (c)  are  generally  colder  than  the  full-sample  means. 
Figures  6,  7  and  8  are  graphical  presentations  of  the  mean  temperatures 
shown  in  parts  (a)  and  (c)  of  Tables  15,  16  and  17. 

Three  factors  should  be  recognized  as  being  relevant  to  the 
July  data  considered.   (i)  In  two  of  the  three  cases,  the  nominal  1% 
extreme  temperature  at  the  forcing  level  T(P  )  did  not  occur,  and  the 
class  was  redefined  by  _a  posteriori  methods  to  comprise  the  set  of  two 
warmest  observations  at  the  forcing  level.   (ii)  The  difference  between 
the  parts  (b)  and  (c)  mean  mandatory  level  profiles  (T  (b)  and  T  (c)), 
at  and  above  300  mb,  never  exceeded  1.0  C.   (iii)  The  dominant  feature 
of  all  of  the  summer  case-studies  was  the  consistent  occurrence  of 
stratospheric-tropospheric  temperature  reversals  noted  in  the 
preceding  paragraph  of  this  section. 
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VI.   CONCLUSIONS 

The  multiple  regression  procedure  presented  for  development  of 
vertically  consistent  model  atmospheres  is  statistically  dependent  upon 
the  input  condition  of  a  known  extreme  temperature  at  a  given  flight 
level.   The  procedure  has  been  shown  to  be  highly  definitive  in  deter- 
mining the  most  probable  vertical  temperature  profile  over  the  stations 
of  interest.   Extremely  warm  stratospheric  temperatures  are  not  asso- 
ciated with  extremely  warm  tropospheric  temperatures,  thus  integrated 
ballistic  density  anomalies  tend  to  be  nearly  compensated.   Likewise, 
extremely  cold  tropospheric  temperatures  tend  to  be  reversed  in  the 
stratosphere.   These  two  points  become  extremely  important  in  the 
determination  of  ballistic  parameters.   As  a  result,  aerothermodynamic 
calculations  can  be  made  for  flight  conditions  involving  changes  in 
altitude  by  using  an  atmospheric  temperature  profile  which  is  consis- 
tent with  the  proposed  MIL-STD-210B  value  at  the  primary  level  of 
interest,  but  is  more  realistic  and  less  severe  at  all  other  levels. 
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